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D . Abstract 

The free energy of a Quark-Gluon Plasma fireball in 

the hadronic medium is calculated in the Ramanathan et 

^ I al. statistical model including the effect of curvature. The 

Li I result with this curvature is found to produce significant 

00 ■ improvement from earlier results in all the parameters we 

"^ ! calculated. The surface tension with this curvature effect 

Ot^^ ' is found to be 0.17 T?, which is two times the earlier value 

QQ ' of surface tension which is 0.078 T^^, and it is nearly close 

O ■ to the lattice value 0.24 T^ . The speed of sound calcu- 

lated with curvature correction is still found to be smaller 
k>( , in comparision with the standard speed of sound in the 



^J^ ■ QGP droplet. 
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1 Introduction 

Quantum chromodynamics predicts [l]a phase transition from a de- 
confined to a confined matter of hadrons during the early stages 
of the universe. Just in the beginning of the universe it was con- 
sidered to be in the form of deconfined matter of quarks and glu- 
ons and subsequent process of coohng leads to transformation into 
hadrons. It is indeed a very complicated phenomena as predicted by 
the heavy-ion collider experiments. So, the formation of quark-gluon 
plasma droplet (fireball) is a very exciting field in the present day 
of heavy ion collider physics [2]. The central assumption in this ap- 
proach treats the QGP-Hadron system as a quasi-static equilibrium 
enabling applicability of equilibrium statistical mechanics to this 
complicated system. This central assumption has been employed by 
many pioneers in this field [2-3]. 

As we know from the extensive literature [4], it is really needed 
to look at the nucleation process produced by the statistical fluctu- 
ations of the critical free energy difference between two phases. The 
model of Csernai-Kapusta et al. [3-5] uses the liquid drop model 
expansion for this nucleation. This model is used by Ramanathan 
et al. [4]. Now it is modified with the effect of curvature term in 
their free energy of the liquid drop. This modified free energy is: 

AF = —R^[P^^,(T,fiB) - PUT.I^b)] 

+ATiR^a + MCR (1) 

The first term represents the volume contribution, the second 
term is the surface contribution where a is the surface tension, and 
the third term is the curvature term. We are interested to see the 
effect of this curvature on the free energy. The critical radius Re 
can be obtained by minimising eqn.(l) with respect to the droplet 
radius -R, which gives two critical radii, 

with smaller radius corresponds to a local mimimum in free energy 
and the larger to a local maximun. We take the local maximum value 
of the radius as suitable case for our fireball system. At the local 
maxima, it has existence of stable solution. Then, we calculate the 
surface tension through AFc = F{Rc, T) — F(0, T) and it is obtained 



i?, = — (1±V(1 + ^^)) (2) 



as: 

Re SttRc 
where ' C ' is the curvature coefficient. To calculate the total rel- 
ativistic density of states for quarks and gluons, it is necessary to 
modify the density of states we obtained in our earlier paper [4]. 
In that paper, we calculated the density of states adapting the pro- 
cedures of the Thomas- Fermi construction of the electronic density 
of states for complex atoms and the Bethe density of states [6] for 
nucleons in the complex nuclei as templates. To modify these to- 
tal density of states, we consider the higher order approximation 
scheme [7] , which can give better results in the calculation of total 
free energy of quarks and gluons. This higher order approximation 
is considered to be the curvature factor in this model. In view of 
this situation, we construct density of states breifly highlighting the 
Ramanathan et al. statistical model in the hadronic medium for 
obtaining free energy of quarks and gluons. In the next section, 
we use total free energy to calculate the interfacial surface tension, 
thermodynamic variables and speed of sound in this model with the 
curvature term. 



2 Detemination of density of states for the QGP 
droplet with curvature term 

To determine the density of states of QGP droplet, the atomic model 
of large atomic number of Thomas and Fermi [6] is modified by Ra- 
manathan et al. [4]. The total electronic density of states defined 
by the Thomas and Fermi in phase space is: 

f Pe{k)dk = [-2mV{k)f^u/37r^ . (4) 

or, 

dV{k) 



Pe{k) = \v{2mfl'llA \-y{ktl^ ■ [-^^J (5) 

This model of electron is replaced by Ramanathan et al.'s model 
of QGP droplet with the corresponding density of states for quark 
and gluon with a suitable QCD induced phenomenological potential 
Vi}^). It is given as: 

^ p^^gdk = [-V,onf{k)fu/37r^ , (6) 

or, 



where u is the volume occupied by the QGP and k is the rel- 
ativistic four-moinentuni in natural units. Vconf{k) could be any 
confining potential for quarks and gluons. This potential plays the 
role of a mean field potential in phase space similar to the mean 
field potential of the Thomas-Fermi scheme, but in a very differ- 
ent context-namely the QGP-Hadron system. The distinction of 
these two model is that Thomas-Fermi model deals about low tem- 
perature feature whereas QGP is high temperature feature which 
subsequently leads to transformation into hadrons. Taking into ac- 
count of all these factors, the density of states was modified with a 
suitable parametrization factor in the dynamics of the QGP fiuid. 
Thus, if we consider that the curvature is essential higher order cor- 
rection factor in free energy expansion, then the density of states is 
further modified by using the technique of Neergaard et al. [7], so 
that it gives more effective results in the free energy. So the modified 
density of states for quarks and gluons is given as: 

P{k) = Pq,g{k) + p{k)cur 

/■ 1 1 
where, p{k)cur = Ci I ds{—- + -—) (8) 

J ri\ 1X2 

here, Ci is function of k/m and depends on the type of field and 
on the boundary conditions. The different values of Ci for different 
particles are given as: 



(10) 

(11) 
which are in the limit of dynamical quark mass. 

3 Mean-field inter-quark potential and the Free 
energy 

To calculate the interacting potential within the system , it is very 
ideal to find the effective mean field potential among the quarks and 
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quark-gluon. This effective potential is obtained through thermal 
hamiltonian for the QGP [8]. 

V^^ik) = {l/2kh,,, g\k)T' . (12) 

where, g{k) for first order QCD running coupling constant, which 
for quarks with three flavors is, [8] 

g\k) = (4/3)(12vr/27){l/ln(l + fcVA')} , (13) 

In the above expression, A = 150MeV is QCD parameter with 
our parametrization factors jg^g which we take 7g = 1/6 and 7g = 8 or 6 
times 7q. This value fits the lattice QCD simulation [9]. This effec- 
tive perturbed potential will have minimum value at each point of 
phase space. That is : 

Vikrmn) = {ig^gN^^T^ k^ /2f'\ (14) 

where iV = (4/3)(127r/27). 

So, this is called as low energy cut off in the model leading to 
finite integrals by avoiding the infra-red divergence. It is of the same 
order of magnitude as A and T. We calculate the free energy for 
quarks and gluons with the above modified density of states. 

F, = TTg., J dkp,{k) ln(l ± e-^V^^^^y^) , (15) 

where Pi{k) is the modified density of states of the particular 
particle i (quarks, gluons, pions etc.) and gi is the degeneracy factor 
(color and particle-antiparticle degeneracy) which is 6 for quarks and 
8 for gluons [7] and 3 for pions. 

The density of state for interfacial surface is calculated through 
the Weyl's model and Ramanathan et al's model. It is obtained as: 

47rR2^ 

Pinter faceyk) — — (16) 

iDTT 

Therefore, the interfacial energy obtained through a scalar Weyl- 
surface in Ramanathan et.al[4, 10-11] with suitable modification to 
take care of the hydrodynamic effects [12] is: 

^interface "7^ V 'J 

where 



7 = v^ X v'(l/7,)2 + (1/7^)2, 



which is the effective rms value of the flow parameter of the 
quarks and gluons respectively. The pion free energy is [10] 

F^ = i'iT/2Tx^)u r edk ln(l - e-V™^+^/^) . (19) 

Jo 

For the quark masses we use the current (dynamic) quark masses 
TTT-o = rrtd = MeV and rris = 150 MeV, just as in reference [10]. 
We can thus compute the total free energy Ffotai as, 

^ total / ^ ^i ~r Winter face ~r ^-k) \^^) 

i 

where i stands for u,d and s quark and gluon. 

4 Interfacial surface tension 

From the above equation (20), we can compute the free-energy 
change with respect to the droplet radius for two promising parametri- 
sation factors. As the value of these two cases give more exciting 
results in the free energy shown in Fig. 1 and Fig.2, we expect that 
the formation of droplet with observable critical droplet radius will 
exhibit the proper transition features in the band of temperature 
150 MeV to 170 MeV which are expected from lattice calculations 
too. 
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Figure 1: Ffotai at 7^ = 67^, 7g = 1/6 for various temperatures. 
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Figure 2: Ftotai at 7g = 87^ , 7^ = 1/6 for various temperatures. 



The Figs. 1 and 2 are the two most promising scenarios of our 
model which exhibit measurable droplet radius of the order of few 
fermi, and show significant barrier amplitudes in free energy which 
account for the nucleation rate of droplet formation, which was 
shown in the earlier paper too [5]. So, we use these promising values 
of the free energy with curvature and look at the changes produced 
by this effect. The set of parameters leading to fig 1 seems more 
realistic because of stability of droplet size at free energy nearly at 
critical radius irrespective of the transition temperature unlike in 
fig. 2. From the values of the critical free-energy at the correspond- 
ing critical fireball radius that can be extracted from the Figs. 1 and 
2, we can compute the surface tension of the fireball using equation 
(3) as listed in Tables 1 and 2. 



T, 


AF, 


lie 


cr 




{MeV) 


(MeV) 


ifm) 


(MeV/fm^) 




150 


921.93 


2.814 


14.798 


0.17 


160 


1059.00 


2.756 


17.949 


0.17 


170 


1200.00 


2.665 


21.528 


0.17 


190 


1486.00 


2.506 


30.083 


0.17 


210 


1765.00 


2.344 


40.586 


0.17 


230 


2030.00 


2.182 


53.365 


0.17 


250 


2272.00 


2.024 


68.540 


0.17 



Table-1 for Surface Tension of QGP droplet at 7^ = 87^,7^ = 1/6. 



T, 


AF, 


Re 


cr 


(J 


(MeV) 


(MeV) 


ifm) 


(MeV/fm^) 




150 


1928.00 


4.416 


9.646 


0.17 


160 


2377.00 


4.488 


18.069 


0.17 


170 


2890.00 


4.555 


21.667 


0.17 


190 


4104.00 


4.659 


30.233 


0.17 


210 


5541.00 


4.713 


40.735 


0.17 


230 


7132.00 


4.703 


53.461 


0.17 


250 


8775.00 


4.618 


68.974 


0.17 



Table-2 for Surface Tension of QGP droplet at 7^ = 67^,7^ = 1/6. 

The constant result of cr/T^ shows that the surface tension is in- 
dependent of transition temperature and the values of the parametri- 
sation factor. 



5 Thermodynamic variables and nature of the 
phase transition 

The thermodynamic properties give the following standard entities 
such as: 

S = -{dF/dT) . (21) 

Cv = T{dS/dT)v . (22) 

Cl = S/Cv (23) 

These entities can be obtained from the total free energy for 
these two promising parametrisation factors jg = 67^ and jg = 
Sjq. The behaviour of entropy, 5" vs temperature, T and the heat 
capacity at constant volume, Cy with temperature, T too indicate 
the real nature of the phase transition of the system. These output 
characteristics are shown in the following Figs. 3 — 8. 
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Figure 3: Variation of S with temperature T at 7^ = 67^ , 7^ = 1/6. 
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Figure 4: Variation of S with temperature T at 7^ = 87, , 7^ = 1/6. 
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The above figs. 3 and 4 indicate that there exists a very weak 
discontinuity in the vicinity of Tc = 160 MeV, in the entropy even 
though the curvature term is included in the total free energy with 
parametrisation factor jg = Qjq and a very mild discontinous in the 
entropy of the parametrisation factor 7^ = 87^. It is happened in 
the first order thermodynamic variable. The discontinuity is just 
of the order of one standard deviation of the entropy variable and 
therefore indeed claim to be a very weak transition . For the second 
order variable Cy shown in figs. 5 and 6, there is absolutely no 
discontinuity in the variable with increase in temperature. 
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Figure 5: Variation of specific lieat Cy with temperature T at 7^ 

67? , Iq = 1/6- 
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Figure 6: Variation of specific lieat Cy with temperature T at 7^ 

87? , Iq = 1/6- 
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Even though the curvature term is accounted for our calculations 
there is a little improvement in the speed of sound. The result is 
rather low in comparision with the lattice calculations of speed of 
sound. The model calculation of sound speed is also independent of 
the transition temperature and model parametrization factor. 



16 



0.4- 



0.3- 



^« 0.2- 
O 



0.1 



0.0- 



***************** — — ♦ * * • 



— I I I I I I I I I I I ^ 
1|Q; 160 180 200 W .340 260 

T(IVleV) 



Figure 7: Variation of speed of sound squared C| with temperature 

T at 7g = 675 , 7g = 1/6. 
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Figure 8: Variation of speed of sound squared C| with temperature 

T at 7g = 875 , 7g = 1/6. 
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6 Conclusion: 

The free energy calculations for the two promising parametrization 
factors are shown in the Figs.[l] and [2]. There is good improvement 
in the free energy with this curvature correction. The energy at Re 
is about two times the energy without the curvature and there is 
increase in the maximum stable droplet size. Due to this significant 
effect in the amplitude of the free energy, we determine the sur- 
face tension, thermodynamic variables and speed of sound in QGP 
for these two factors. There is significant increase in surface ten- 
sion with the critical temperature in this modified free energy. This 
shows good improvement over the earlier results of Ramanathan et 
al. without curvature. This leading order term of curvature has 
good impact to find the surface tension. The calculated value of 
surface tension is found to be 0.17 T^f , which is two times the earlier 
results that is 0.078 T^ [5] and this present value is nearly close to 
higher value of lattice results that is < 0.24T^''[7,13]. The result is 
still consistent with earlier ones in terms of its constancy through 
out the temperature and parametrization factors. It is perfect con- 
formity with the latest QCD simulations [9]. The model still shows 
weakly first order phase transition at the temperature in the range 
(160 ± 5) MeV with this effect of curvature in the free energy 
and it is expected with the current feature of QGP-Hadron phase 
transition too [1]. Thus, this transition is a mild characteristic in 
nature and discontinuity is found in the first order thermodynamic 
variable namely the entropy S as shown in Figs. 3 and 4. Lastly, 
the speed of sound with curvature correction is compared with the 
earlier result without curvature. The improvement in the speed of 
sound is very small . We obtained this value earlier nearly close to 
0.09 T jTc- Presently it is found to be 0.1 T jT^ but still the result 
is showing low value with the lattice QCD as before[13]. This low 
result is expected with the recent model calculations [14] and lattice 
simulations that include dynamical quarks [15], but as should be 
expected, it is at variance with pure gauge lattice results[16]. 
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